Complex Phase Behavior of Repulsive Step Potential Systems in Three Dimensions 
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The comprehensive computer simulation study of the phase diagram of the repulsive step potential 
system in three dimensions is represented. We show that the system with a simple purely repulsive 
isotropic potential demonstrates a number of unusual features. The maxima and minima on the 
melting curve are found for some regions of potential parameters. It is shown that the phase 
diagram in p — T plane includes two isostructural crystalline parts separated by the disordered 
phase which is amorphous at low enough temperatures. Phase diagram in the (P-T) plane shows 
that the transition to the amorphous state occurs approximately along the extrapolated spinodals. 
Structural FCC-BCC phase transition is found at high densities. 

PACS numbers: 61.20.Gy, 61.20.Ne, 64.60.Kw 
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After the pioneering work by Hemmer and Stell 
where the soft core potential with an attractive interac- 
tion at large distances was first proposed for the qualita- 
tive explanation of the isostructural phase transitions in 
materials such as Ce or Cs, a lot of attention was paid 
to the investigation of the properties of the systems with 
the potentials that have a region of negative curvature in 
their repulsive core [2j, y, U, la, la Lfl, la, M, LL0, llll, iTJ, 1 14 
These studies are particularly relevant 
because the interatomic potentials of some pure metal- 
lic systems, metallic mixtures, electrolyte and colloidal 
systems can be approximated by such type of potentials. 
The simplest form of the negative curvature potential is 
the repulsive step potential which consists of a hard core 
plus a finite repulsive shoulder of a larger radius. It is 
well known that systems of particles interacting through 
the pair potentials which have repulsive shoulders can 
possess a rich variety of phase transitions and thermo- 
dynamic anomalies, including water-like anomalies 
liquid- liquid phase transitions [a, la Ml and isostructural 
transitions in the solid region [16|, |l7J, ll8[ . The features of 
greatest interest in the experimental phase diagrams of 
Ce and Cs are the isostructural solid-solid phase transi- 
tions and the melting temperature minima and maxima. 
The authors of Ref. 2, 17 1 succeeded in modelling these 
features of the phase diagram in two dimensions. 



In this paper we present results, based on computer 
simulations in three dimensions, of the phase behavior 
of the system with the purely repulsive step potential 
which reveal a surprisingly complex phase behavior of 
the system as a function of the potential parameters and 
show new features not found in previous studies. 



The repulsive step potential has the form: 
$(r) = 




(1) 



and is characterized by the hard core diameter d, the 
width of the repulsive step a and its height e. In the 
low temperature T « e and high temperature T >> e 
limits the system behaves like a simple hard sphere sys- 
tems with hard-sphere diameter a or d correspondingly. 
Schematic phase diagram may be understood in the fol- 
lowing way: the melting line of the hard sphere system 
has the form [2(| P = cT/a' 3 , where c « 12 and a' is a 
hard sphere diameter. At low and high temperature lim- 
its one has two straight melting lines which correspond to 
hard spheres with diameters a and d. At higher temper- 
atures and densities the line corresponding to a should 
bend to the line corresponding to the hard spheres with 
the diameter d. These lines merge at T « e and, de- 
pending on the magnitude of the ratio s = a/d one can 
expect more or less pronounced anomalies in this region. 
For high values of s one can expect well pronounced max- 
imum on the resulting melting curve which should disap- 
pear in the limit s — > 1 19). 

In our simulations we have used a smoothed version of 
the repulsive step potential (JTJ> , which has the form: 



d>(r) = - 
r 



tanh(fc n (r-CT s ))) (2) 



where n = 14, ko = 10. We have considered the fol- 
lowing values of a s : a s = 1.55,1.35,1.15. Further in 
this paper we use the dimensionless quantities: f = r/d, 
P = Pd 3 /e, V = V/Nd 3 = 1/p, f = k B T/e, omitting 
the tilde marks. 
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FIG. 1: A repulsive step potential consisting of a hard core 
plus a finite shoulder (dashed line) (e = 1, a = 1.5) along 
with the continuous version of the potential ([2]) used in the 
simulations (s = 1, a s = 1.55). 
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FIG. 3: RDF (upper figure) and MSD for different densities 
at T = 0.02. 
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FIG. 2: The family of pressure isotherms for the purely re- 
pulsive step potential ((2} for a 3 — 1.55. 



In Fig. Q] the repulsive step potential is shown along 
with its smooth version which was used in our Monte- 
Carlo and molecular dynamics simulations. 

In order to investigate phase diagrams of the system 
with the potential © we studied isotherms, isochores, ra- 
dial distribution functions (RDF) , mean-square displace- 
ments (MSD), positional order parameter (POP) in the 
framework of the standard molecular dynamics in the 
NVT ensemble with periodic boundary conditions. Tem- 
perature is fixed by rescaling the velocities of the parti- 
cles whenever necessary [2l|. Pressure is calculated by 
a direct evaluation in terms of intcrparticle forces. To 
obtain RDF, MSD and POP, we started from the lat- 
tice structure (FCC or BCC), equilibrated for 5 x 10 6 
cycles and than measured for 3 x 10 6 cycles with a time 
step 5t = 5 x 1CP 5 , or from disordered phase obtained by 
starting from the crystal one at high temperature and low 
density. In some cases the system demonstrates glassy 
behavior therefore equilibration run had to be so long. 

We also calculated Helmholtz free energy of the sys- 



tem using the thermodynamic integration in the frame- 
work of NVT ensemble Monte Carlo simulations 2l| . For 
calculating the free energy of a liquid phase we used 
integration from dilute gas limit [21|. Solid phase free 
energy was calculated by coupling with Einstein crystal 
(Frenkel-Ladd method) [2l],[22|]. We did 10 6 equilibration 
cycles and then 50000 production cycles. 

In order to locate roughly phase boundaries for crystal 
phases we analyzed the behavior of RDF, MSD and POP 
in the vicinity of van der Waals loops on the isotherms 
using MD simulations. In Fig. [5] the family of pressure 
isotherms for the purely repulsive step potential @ with 
a s = 1.55 for different values of T is shown. One can see 
that for T = 0.02, 0.05, 0.1 there are 3 sharp bends on the 
pressure isotherms and only one sharp bend for T = 0.2. 
The nature of the transitions may be understood from the 
behavior of RDF and MSD shown in Fig. |3] for T = 0.02. 
It can be seen that for p = 0.15 (before the first bend 
on the corresponding isotherm) RDF and MSD have a 
liquid-like character. For p — 0.2 RDF and MSD have 
solid-like forms. The most interesting behavior corre- 
sponds to the densities in the range between second and 
third bends on the isotherms. In this case (see lines cor- 
responding to p = 0.5) RDF is liquid-like, but MSD has 
solid-like form. This means that in the range of densi- 
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FIG. 4: Phase diagram of the system of particles interacting 
through the potential (2) with o s — 1.55 in p — T and P — T 
planes. 



ties between second and third bends on the isotherms at 
low enough temperatures there exists an amorphous solid 
phase, which transforms to the crystalline solid after the 
third transition (see RDF and MSD for p — 0.8). In this 
density region the crystalline solid is unstable against the 
amorphization at all temperatures. The lowest tempera- 
ture at which the simulations were done was T = 0.00001. 
However, it may be shown that at high temperatures the 
system is a liquid (RDF and MSD are liquid- like) . The 
location of the liquid-glass transition is a subject of the 
separate publication. This approach gives the possibility 
to obtain the approximate limits of stability (spinodals) 
of solid and liquid phases. 

In Fig. QJa) we show the phase diagram of the sys- 
tem with a s — 1-55 as obtained using the double tangent 
construction. Open circles determine the boundaries of 
liquid and solid (FCC) phases. Squares correspond to 
the structural FCC-BCC transition. As it was discussed 
above there are two FCC phases separated by the rather 
wide density range of the disordered phase which is amor- 
phous at low temperatures. To our knowledge this is a 
new type of phase behavior which was not found in pre- 
vious studies. 

In Fig. |IJb) the phase diagram of the system with 



<7 S = 1.55 is represented in (T — P) plane. The cir- 
cles correspond to the melting line obtained from the 
double tangent construction, and open triangles denote 
the spinodals obtained from MD simulations. As one 
can expect the spinodals lie higher in temperature than 
the thermodynamic transition. At higher densities there 
is a phase transition between FCC and BCC phases 
(squares). The line of the FCC-BCC transition was ob- 
tained using the double tangent construction. MD simu- 
lations (RDF and POP behavior) confirm these results. 
The thermodynamic melting line (circles) consists of two 
parts which should merge at low temperatures. Unfortu- 
nately, we could not calculate the phase diagram between 
these parts because of problems with calculation of the 
Hclmholtz free energy of the disordered state. 

In the insert in Fig. |U(b) the stability line of the crys- 
talline solid phase is shown by the squares. The open 
circles correspond to the stability limits of the disor- 
dered phases obtained from the pressure isotherms (see 
Fig. [2|). One can see that the transitions to the amor- 
phous phases occurs along the extrapolations of the melt- 
ing lines to the supercooled region, and the " melting" in 
this region is nothing but the pressure induced amor- 
phization 



14, 



\m \M I2J, 125J. It should be noted that 
due to the fact that only spinodals are shown in the in- 
sert in Fig. 0Kb) , we can not determine whether pres- 
sure induced amorphization is conventional "two-phase 
melting" or mechanical instability [3, 15, [23[. In or- 



der to explain the phase diagram shown in the insert in 
Fig. [U(b) let us consider the behavior of the system under 
compression at temperatures lower than the melting line 
minimum. Taking into account Fig. [21 it can be seen that 
upon increasing pressure one intersects the left hand part 
of curve 1 corresponding to the stability line of the liquid 
phase. At this line liquid-low-density FCC solid transi- 
tion occurs. At the continuation of line 2 the transition 
from a low-density FCC solid to amorphous phase takes 
place. Line 4 corresponds to the limit of stability of the 
amorphous phase. At this line amorphous phase crystal- 
lizes and becomes a high-density FCC phase. Upon de- 
compression from the high-density FCC phase one first 
meets the continuation of line 3 where the amorphiza- 
tion occurs. Right hand part of line 1 corresponds to the 
limit of stability of the amorphous phase, where amor- 
phous phase crystallizes to the low-density FCC phase. 
Upon further decompression the FCC phase melts at the 
left hand part of line 2. 

It should be emphasized that we can not reach true 
thermodynamic equilibrium state in the density range 
where we find the amorphous phase. In principle, the 
energy of the amorphous phase may be higher than the 
energies of FCCi and FCC2 solids, and in this case the 
amorphous phase may exist only because of the very slow 
nonergodic kinetics of the system. If the system could be 
equilibrated one would find the isostructural FCC1-FCC2 
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FIG. 5: Phase diagram of the systems of particles interacting 
through the potential © with a s = 1.15, 1.35, 1.55 in T - P 
plane. Filled circles correspond to the melting line of the soft 
sphere potential "I>(r) = 1/r 14 without repulsive shoulder. 



phase transition (see, for example, [181]'). 

In Fig. [5] the evolution of the melting line obtained 
by equating the Helmholtz free energies is shown as a 
function of the size of the repulsive step er s . It can be 
seen that the system behaves in qualitative accordance 
with the described above model of the mixture of the 
hard spheres with two diameters (l9j . In particular, the 
melting curve maximum disappears for a s = 1.15. We 
are going to discuss this phase diagram in more details in 
a separate publication. For the comparison we also show 
the melting line of the soft sphere potential = 1/r 14 
without repulsive shoulder (filled circles). 

In conclusion, we performed the comprehensive com- 
puter simulation study of the phase diagram of the sys- 
tem described by the simple purely repulsive step po- 
tential @ in three dimensions and found a surprisingly 
complex phase behavior. The maxima and minima on 
the melting curve are found for some regions of potential 
parameters. It is shown that the evolution of the phase 
diagram with change of the parameters may be qualita- 
tively explained in the framework of the simple model in 
which the system is considered as a mixture of two types 
of hard spheres. For the first time it is shown that the 
phase diagram in p — T plane includes two FCC phases 
with different densities separated by the disordered phase 
which is amorphous at low enough temperatures. We 
show that the transition to the amorphous state occurs 
approximately along the extrapolated spinodals in the 
(P-T) plane. Structural FCC-BCC phase transition is 
found at high densities. 

We thank V. V. Brazhkin for stimulating discussions. 



[1] 
[2] 

[3 
[4 
[5 
[6 
LX 
[8 
[9 

[io; 
in 

[12 
[13 
[14 

[is; 

[16 

[17; 
[is; 

[19 
[20 

[21 

[22 

[23 
[24 

[25] 



P. C. Hemmer and G. Stell, Phys. Rev. Lett. 24, 
1284(1970); G. Stell and P. C. Hemmer, J. Chem. Phys. 
56, 4274 (1972). 

V. V. Brazhkin. S. V. Buldyrev, V. N. Ryzhov, and H. E. 
Stanley [eds], New Kinds of Phase Transitions: Trans- 
formations in Disordered Substances [Proc. NATO Ad- 
vanced Research Workshop, Volga River] (Kluwer, Dor- 
drecht, 2002). 

E.Velasko, L. Mederos, G. Navascues, P. C. Hemmer, and 

G. Stell, Phys. Rev. Lett. 85, 122 (2000). 

P. C. Hemmer, E.Velasko, L. Mederos, G. Navascues, and 

G. Stell, J. Chem. Phys. 114, 2268 (2001). 

V. N. Ryzhov and S. M. Stishov, Zh. Eksp. Teor. Fiz. 

122, 820 (2002)[JETP 95, 710 (2002)]. 

V. N. Ryzhov and S. M. Stishov, Phys. Rev. E 67, 

010201(R) (2003). 

Yu. D. Fomin, V. N. Ryzhov, and E. E. Tareyeva, Phys. 
Rev. E 74, 041201 (2006). 

M. R. Sadr-Lahijany, A. Scala, S. V. Buldyrev and H. E. 

Stanley, Phys. Rev. Lett. 81, 4895 (1998). 

M. R. Sadr-Lahijany, A. Scala, S. V. Buldyrev and H. E. 

Stanley, Phys. Rev. E 60, 6714 (1999). 

P. Kumar, S. V. Buldyrev, F. Sciortino, E. Zaccarelli, 

and H. E. Stanley, Phys. Rev. E 72, 021501 (2005). 

L. Xu, S. V. Buldyrev, C. A. Angell, and H. E. Stanley, 

Phys. Rev. E 74, 031108 (2006). 

E. A. Jagla, J. Chem. Phys. Ill, 8980 (1999). 

E. A. Jagla, Phys. Rev. E 63, 061501 (2001). 

S. Bustingorry and E. A. Jagla, Phys. Rev. B 69, 064110 

(2004) . 

S. Bustingorry and E. A. Jagla, Phys. Rev. B 71, 224119 

(2005) . 

D. A. Young and B. J. Alder, Phys. Rev. Lett. 38, 1213 
(1977). 

D. A. Young and B. J. Alder, J. Chem. Phys. 70, 473 
(1979). 

P. Bolhuis and D. Frenkel, J. Phys.: Condens. Matter 9, 
381 (1997). 

S. M. Stishov, Phil. Mag. B 82, 1287 (2002). 

W. G. Hoover, S. G. Gray, and K. W. Johnson, J. Chem. 

Phys. 55, 1128 (1971). 

Daan Frenkel and Berend Smit, Understanding molecu- 
lar simulation (From Algorithms to Applications), 2nd 
Edition (Academic Press), 2002. 

D. Frenkel and A. J. C. Ladd, J. Chem. Phys. 81, 107 
(1984). 

O. Mishima, Nature (London) 384, 546 (1996). 

A. G. Lyapin and V. V. Brazhkin, Phys. Rev. B 54, 12 

036 (1996). 

E. L. Gromnitskaya, O. V. Stalgorova, V. V. Brazhkin, 
and A. G. Lyapin, Phys. Rev. B 64, 094205 (2001). 



